Primordial Non-Gaussianity in the Forest: 3D Bispectrum of Ly-a Flux Spectra 

Along Multiple Lines of Sight 
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We investigate the possibility of constraining primordial non-Gaussianity using the 3D bispectrum 
of Ly-a forest. The strength of the quadratic non-Gaussian correction to an otherwise Gaussian 
primordial gravitational field is assumed to be dictated by a single parameter / NL . We present 
the first prediction for bounds on / NL using Ly-a flux spectra along multiple lines of sight. The 
3D Ly-a transmitted flux field is modeled as a biased tracer of the underlying matter distribution 
sampled along ID skewers corresponding to quasars sight lines. The precision to which / NL can 
be constrained depends on the survey volume, pixel noise and aliasing noise (arising from discrete 
sampling of the density field). We consider various combinations of these factors to predict bounds 
on / NL . We find that in an idealized situation of full sky survey and negligible Poisson noise one may 
constrain / NL ~ 23 in the equilateral limit. Assuming a Ly-a survey covering large parts of the sky 
(kmin — 8 x 10~ 4 Mpc _1 ) and with a quasar density of n = 5 x 10 -3 Mpc -2 it is possible to constrain 
/ NL ~ 100 for equilateral configurations. The possibility of measuring / NL at a precision comparable 
to LSS studies maybe useful for joint constraining of inflationary scenarios using different data sets. 

PACS numbers: 98.80. Cq, 98.62.Ra 



Introduction: The widely popular paradigm of slow-roll 
inflation driven by a single canonical scalar field gener- 
ates adiabatic perturbations which are largely Gaussian 
in nature and leads to a nearly scale invariant power spec- 
trum [![. Several theoretical predictions however point 
towards mild to severe departure from Gaussianity |2|. 
Measuring the degree of non-Gaussianity is hence crucial 
towards discriminating between various inflationary sce- 
narios thereby enhancing our understanding of the very 
early Universe. It is assumed that on sub-horizon scales 
the primordial gravitational potential $P rim is related to 
a Gaussian random field 4>g through a non linear rela- 
tion of the form $P rim = 4> G + (<j) G 2 - (0 G 2 )). De- 
partures from Gaussianity may hence be quantified using 
the parameter / NL . We assume that / NL is scale indepen- 
dent - a reasonable prediction for most inflationary mod- 
els where non-Gaussianity is generated on super-horizon 
scales. The value of f„, obtained from slow-roll infla- 
tion turns out to be very small (O (10~ 3 — 10~ 2 )) [l|. 
This implies that any detection of large / NL shall rule out 
all canonical single field slow-roll inflation models. The 
mean value of / NL (26 ± 140 for equilateral and 32 ± 21 
in local limit for la CL) obtained from WMAP data Q 
seems to indicate large non-Gaussianity. Although low 
SNR in these results indicate that we are yet to detect 
the primordial non-Gaussianity, it is expected that data 
from the Planck satellite 1 shall give a much tighter con- 
straint on / and the error is expected to come down to 



±5 in the local limit. Other than CMBR obser- 



vations, a measurement of the bispectrum or the three 
point correlation function of the galaxy distribution [3] 
is a standard alternative method to constrain primor- 
dial non-Gaussianity. These probes however only provide 
weak bounds on / NL as compared to the CMBR observa- 
tions (SDSS 2 for example can measure |/ NL | ~ 10 3 — 10 4 ). 

In the post reionization epoch, small fluctuations of the 
neutral hydrogen (HI) density field in a predominantly 
ionized IGM leads to a series of distinct absorption fea- 
tures - the Ly-a forest Q in the spectra of background 
quasars. The Ly-a forest is a well established power- 
ful probe of cosmology 0, 0j- Traditional Ly-a stud- 
ies have considered the power spectrum or bispectrum 
of the one dimensional transmitted flux field [7| corre- 
sponding to the quasar line of sight. This approach is 
reasonable when the angular density of quasars on the 
sky is low. The new generation of quasar surveys (the 
ongoing BOSS 3 and future BigBOSS 4 ) however promise 
to achieve a very high quasar density and cover large 
fractions of the sky. This has led to the possibility of 
measuring the 3D Ly-a power spectrum along multiple 
lines of sight [8| . We note that the first hydro simulation 
of the Ly-a forest in non-Gaussian scenarios is presented 
in an earlier work [9|. 

In this letter we investigate the possibility of constrain- 
ing / NL using the 3D Ly-a forest bispectrum. Similar to 
the power spectrum studies, the Ly-a flux distribution is 
assumed to be a biased tracer of the underlying matter 
field sampled along discrete sight lines. We explore the 
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range of observational parameters for the constraints on 
/ NL from the 3D analysis to be competitive with CMBR 
and LSS (large scale structure) studies. 

Formalism: The post-reionization matter density field 
A k in Fourier space is related to the primordial grav- 
itational potential on sub- horizon scales as Ak(z) = 
M(k,z)$l lim . The function M(k,z) is given by 

M(k, z) = — | Q T gi D+ (z), where T(k) denotes the mat- 
ter transfer function and D+(z) is the growing mode of 
density fluctuations. We have used the BBKS trans- 
fer function [Io[ and cosmological parameters ob- 
tained from a MCMC 5 analysis on the WMAP7 data- 
sets. The power spectrum and bispectrum of the density 
field are defined as (A kl Ak 2 ) — <5o(ki + k 2 )P(fci) and 
(A kl A k2 A k3 ) = 8 D (ki +k 2 +ka)-B(Ai,*2,*3). Clearly, 
the linear power spectrum of the density field is given by 
P L (fc) = A4(k, z) 2 P|) rim where P<^ rim denotes the primor- 
dial power spectrum of the gravitational potential such 
that P| rim = P 0G + 0(f*J. The powerspectrum P 0G of 
the Gaussian field G shall be assumed to be featureless 
and scale invariant. We note that a wide class of inflation 
models, including the simplest one comprising of a single 
inflaton field in a quadratic potential, introduces pertur- 
bations that are almost Gaussian and exhibit a power 
spectrum that is nearly scale invariant. It follows that in 
the realm of linear perturbation theory the bispectrum of 
the matter field arising from primordial non-Gaussianity 
is given by P L i 2 3 = M(h)M (k 2 )M(k 3 )B^ Gl23 where 
we use the notation 123 = {k\,k 2 ,k 3 ) and B^q is given 
by 



B 



<#> G 123 



2/ NL 



[P0 G (fci)P 0G (fc 2 ) + eye] + 0(/ 3 L ). (1) 



Apart from the contribution to the bispectrum from 
primordial fluctuations, non-linear structure formation 
caused by gravitational instability leads to mode coupling 
and thereby induce additional non-Gaussianity. This is 
especially relevant when we use low redshift tracers to 
implicitly measure n— point functions of the matter den- 
sity field. Using the second order perturbation theory 
the additional contribution to the matter bispectrum is 



nNL 
13 123 



2P 2 (k 1 ,k 2 )P(fc 1 )P(fc 2 ) + cyc. 



(2) 



where the second order correction is obtained from per- 
turbation theory and we have adopted the form of P 2 
from [4]. Filially, the total matter bispectrum is a sum 
of the contributions to non-Gaussianity arising from the 
intrinsic primordial fluctuations and that generated by 
non-linear evolution of an otherwise Gaussian field. Thus 
we have Pi 2 3 = B\ 2 3 + Bf^, where we have ignored the 
possible contribution from the primordial trispectrum. 



We use Pi 2 3 to obtain the 3D bispectrum of the Ly-a 
forest. 

The Ly-a forest spectra are associated with gas dis- 
tribution in voids or slightly overdense regions. Not- 
ing that the astrophysical structures associated with the 
spectra are only mildly non-linear, the transmitted flux 
T through the Ly-a forest may be modeled by assuming 
that the gas traces the underlying dark matter distri- 
bution 0] except on small scales where pressure plays 
an important role. Further, it is believed that photo- 
ionization equilibrium that maintains the neutral frac- 
tion also leads to a power law temperature-density re- 
lation (l2| . The fluctuating Gunn- Peterson approxima- 
tion (FGPA) [l3| incorporates these assumptions to relate 
the transmitted flux T to the dark matter overdensity 
S as T = T exp [-A(l + £) 2 -0-7(7-i)] ; where j is thc 

mean transmitted flux and (7 — 1) is the slope of the 
temperature-density relation [12j. We note that 7 im- 
prints the reionization history of the Universe. The red- 
shift dependent quantity A [l4| depends on a number of 
parameters, like the IGM temperature, photo-ionization 
rate and cosmological parameters Q- On large scales 
it is reasonable to believe that the fluctuations in the 
transmitted flux 8jr = {TjT— l) may be expanded as 
Sjr = b\5 + \b 2 5 2 where, it is assumed that the Ly-a for- 
est spectrum has been smoothed over some suitably large 
length scale. This relation allows analytic computation 
of statistical properties of 8p. We note that corrections 
to this on small scales come from peculiar velocities, an 
effect we have not incorporated in our analysis for sim- 
plicity. At our fiducial redshift z — 2.5, we adopt an 
approximate (T, 7, A) = (0.8, 1.5, 0.16) from the numer- 
ical simulations of Ly-a forest [15| and theoretical pre- 
dictions 16[. We note however that these numbers are 



largely uncertain owing to inadequate modeling of the 
IGM. The bias, b\ for example has a sensitive redshift 
dependence and may depend on the smoothing scale of 
the Ly-a spectra. Using the local bias model, the power 
spectrum Pr(fc) and bispectrum Bjr of Ly-a forest flux 
fluctuations Sjr are given by 

Pr(fc) = bjP(k) 

Bru3 = b\B l2Z + b\b 2 [P{k x )P{k2) + eye] (3) 



The bispectrum of Ly-a flux is hence completely modeled 
using three parameters (/ NL , b±, b 2 ). We shall now set up 
the Fisher matrix for constraining / NL using the Ly-a 
bispectrum. Following the formulation described in [3] 
we define the bispectrum estimator as 
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> 123 = T7— \ d 3 qi / ^ 3 q2 / d 3 q3<5z)(qi23) 

V123 Jk! Jk 2 Jk 3 

x A^(fc 1 )A^(fc 2 )A^(fc 3 ) (4) 

Here qi 2 3 = qi + q 2 + q3, Vf = (27r) 3 /U where V is the 
survey volume, Ui 23 = J*^ d 3 qi d 3 q 2 J fca d 3 q 3 <5z)(qi23) 
and the integrals are performed over the qi— intervals 
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(h - $£,ki + 4p). The quantities A^Qm) denotes the 
'observed' Ly-a flux fluctuations in Fourier space. The 
observed quantity Sjr(r) is given by the continuous field 
Sjr(r) sampled along skewers corresponding to line of 
sight to bright quasars. We therefore have <5£-(r) = 
(5jr(r) x p(r), where the sampling window function p(r) 

TV" ^" w a S^(r ± -r ±a ) where N . s a 



is defined as p(r) 



El 



normalization such that J dVp(r) = 1. The summation 
extends up to Nq, the total number of quasar skewers in 
the field which are assumed to be distributed with sky 
locations r± a . The weights w a introduced in p(r) are in 
general related to the pixel noise and can be chosen with a 
posteriori criterion of minimizing the variance. In Fourier 
space, we then have A^k) = /5(k)<g>Ajr(k) + Ajr n oisc(k), 
where p is the Fourier transform of p, and A_F noise (k) de- 
notes a possible noise term. 

If the bispectrum covariance matrix is diagonal which 
implies that no correlation exists between different trian- 
gle shapes, the simple variance of the estimator Bjr can 

~ 2 -2 

be calculated as ABjr = {B T ) — (Bjr) 2 . This is given 
at the lowest order by 

Vf 



AS, 



^123 



(fc 1 )PJ ot (fc 2 )PJ ot (fc 3 ) 



(5) 



where s = 6, 1 for equilateral and scalene triangles re- 
spectively and Pj ot (k) is the total power spectrum of 
Ly-a flux given by 

(6) 



P><*(k) = P^k) 

>1D/ 



P 



ID/ 



N 



(fen)JV 

The quantity PjP(k\\) is the usual ID flux power spec- 
trum Q corresponding to individual spectra given by 
pi. D (fc||) = (2ir)- 2 f d 2 k ± Pjr(k) and P w denotes the 
power spectrum of the window function. The quantity 
Njr denotes the effective noise power spectra for the Ly- 
a observations. The term PjP(k\\)Pw referred to as the 
'aliasing' term, is similar to the shot noise in galaxy sur- 
veys and quantifies the discreteness of the ID Ly-a skew- 
ers. It has been shown that an uniform weighing scheme 
suffices when most of the spectra are measured with a 
sufficiently high SNR This gives Pw = where n 
is the 2D density of quasars (n = Nq/A, where A is 
the area of the observed field of view). We assume that 
the variance c^at °f the pixel noise contribution to Sjr is 
the same across all the quasar spectra whereby we have 
Njr = a^r N /n for its noise power spectrum. In arriving 
at equation (|6]) we have ignored the effect of quasar clus- 
tering. In reality, the distribution of quasars is expected 
to exhibit clustering (l7| . However, for the quasar sur- 
veys under consideration, the Poisson noise dominates 
over the clustering and the latter may be ignored. 

The Fisher matrix for a set of parameters pi is con- 
structed as 

k 2 



F, 



ki—k mi 



fci 
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k2=k m -, 
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1 dB 



T123 



dB 



T123 



■ k 3 =k„ 



AB F 2 d P* d Pj 



where k m i n = max(k m i n , \ k± — A^l) and the summa- 
tions are performed using 5k = k m i n . Assuming the like- 
lihood function for pi to be a Gaussian the errors in p t 
is given by the Cramer-Rao bound of = F^ 1 . We have 
used this to investigate the power of a Ly-a survey to 
constrain / NL . 

Results: We consider quasars in the range z — 2 to 3 
since the peak in redshift distribution of quasars occur in 
this range (lij . We note that for a given quasar at red- 
shift z = zq, the proximity effect will not allow the spec- 
trum to be measured in the region 10, 000 km s -1 blue- 
wards of the Ly-a emission and only the region which is 
at least 1,000 kms -1 red-ward of the quasar's Ly-/3 and 
O-VI lines are considered to avoid the possible confusion 
with these lines. We have chosen z — 2.5 as out fidu- 
cial redshift for the subsequent analysis. We note here 
that all the parameters involved in the modeling the Ly-a 
forest, have direct or indirect redshift dependence. 

A Ly-a forest survey towards measurement of power 
spectrum or bispectrum is characterized by the survey 
volume, pixel noise in the spectra and the number den- 
sity of the quasar skewers. The constraining power of 
the survey shall depend directly on the choice of these 
parameters. In the cosmic variance limit the minimum 
/ NL that can be measured depends on the number of 
Fourier modes in the survey volume V given by Np- — 
4:TT/3k 2 nax V/(2Tr) 3 . Clearly the minimum detectable / NL 
is a function of k max and k m i n . The noise power spec- 
trum Njr is given by Njr = J" 2 [S*/iV]^(Aa;/lMpc) 
where [S/N]a x is the signal to noise ratio for a spec- 
trum smoothed to a resolution Ax. We quote [S/N] here 
for 1 A pixels. The main source of noise to the 3D power 
spectrum comes from the aliasing noise term and one re- 
quires a very high density of quasars in the field of view 
for this term to be sub-dominant. 
The bispectrum SNR depends on the triangle config- 
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analysis for various combinations of (k m in, n, S/N). 



urations considered to evaluate it. In this work we 
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FIG. 1: The 68.3%, 95.4% and 99.8% likelihood confidence 
contours for the parameters (/ NL , bi). Shown in the figure are 
the values (k m i n ,fi, S/N) used to compute the Fisher matrix. 



have used the simplest equilateral configurations char- 
acterized by just a single Fourier mode. This over esti- 
mates the noise by at least a factor of ~ 2.45 as com- 
pared to the case with arbitrary triangles. In the equi- 
lateral limit the 3D Ly-a bispectrum can be written as 



B F {k) = P{k) 
a 2 = &b\F 2 



M{k) 



3bfb 2 . 



a 2 where a\ = 66^/ NL /c 2 and 

Only two parameters are sufficient 
to model the bispectrum instead of three parameters 
(/ NL , bi, b 2 ) for the general case. We use the fiducial val- 
ues (/ NL , b\, b 2 ) = (0, —0.15, —0.075) and choose / NL and 
bi to be the free parameters for the Fisher analysis. We 
recall that in our modeling of the Ly-a forest we used the 
parameters (J 7 , A, 7). The parameter J 7 does not appear 
in Sjr and there is degeneracy between the parameters A 
and 7 which only appears as a product in b±. Changing 
b\ hence amounts to changing either or both A and 7. 

We assume that the likelihood function is a bivariate 
Gaussian which yields the confidence ellipses shown in 
figure ([T]). The tilt of the error ellipses indicate corre- 
lation between the parameters. We quantify this using 

1 For the 
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the correlation coefficient r = F : 

range of parameters chosen we find that this is roughly 
constant r ~ —0.7. 

In the ideal situation of full sky coverage and negligible 
Poisson noise we find that A/ NL ~ 23 in the equilateral 
limit. We tabulate our results for varying sky coverage 
(^min = ^/(^t 1 ") 3 ), Poisson noise (~ 1/n) and pixel noise 
{S/N) in TableQ] As expected we have tighter constraints 
on (/ NL ,&i) with increasing survey volume, n and S/N. 
The values of the survey parameters chosen are reason- 
able and achievable by future Ly-a surveys. Exploiting 
the entire sky coverage of SDSS we find that one can ob- 



tain a bound on / N 
for a survey with n = 5 x 10" 
are measured at 5a level. 

Our analysis has largely focussed on the equilateral 
configuration. However we find that the Cramer-Rao 
bound for /jvl in the squeezed limit (£3 << rain(k\, k 2 )) 
turns out to be ~ 40 — 100 for the cases we have consid- 
ered. The case of arbitrary triangular configuration is 
to be addressed in our future work [20(. However, our 
preliminary estimates show us that we may constrain 
Jnl ~ 1 in an ideal environment. For example using 
S/N ~ 5, n ~ 10" 3 Mpc" 3 and k mm ~ 10 -3 Mpc _x , we 
have A/jvl ~ 5 in the case of arbitrary triangles, which 
is competitive with CMBR and LSS studies. 

Our discussion so far has bypassed the issue of red- 
shift space distortion arising from peculiar motion. The 
bispectrum in redshift space shall depend not only on 
the shape of the triangle but also on its orientation 
with the line of sight. This shall introduce two more 
angular parameters in the analysis. The multipole ex- 
pansion of bispectrum in redshift space provides a way 
to break the degeneracy between bias and cosmological 
growth parameter / ~ ^m 6 - We plan to take this up 
in a future work [20j . However the spherically averaged 
bispectrum is related to its real space counterpart as 
B%r = (1 + 4/3 + ±/3 2 ) where = f/b t 19]. This 



is an enhancement on all scales 



To conclude, we emphasize that it is possible to put 
stringent bounds on primordial non-Gaussianity from the 
measured 3D bispectrum of the Ly-a forest along multi- 
ple lines of sight and thereby constrain various inflation- 
ary scenarios. Our analytic predictions indicate that such 
studies with future Ly-a surveys may be useful while per- 
forming a joint analysis using other data sets like CMBR 
or LSS. 
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